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O ■ ABSTRACT 

(N 

»— 5 \ We compute locations of sonic points and standing shock waves in a thin, ax- 

CH \ isymmetric, adiabatic flow around a Schwarzschild black hole. We use completely 

analytical method to achieve our goal. Our results are compared with those ob- 
tained numerically and a good agreement is seen. Our results positively prove 
the existence of shocks in centrifugal pressure dominated flows. We indicate how 
our results could be used to obtain spectral properties and frequencies of shock 
oscillations which may be directly related to the quasi-periodic oscillations of 



hard X-rays. 
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1. Introduction 

>' 
"^ ■ In recent years, study of standing and oscillating shocks in accretion flows has become 

very important since it is recognized that the spectral states of black holes as well as Quasi- 
Periodic Oscillations (QPOs) observed in light curves of black hole candidates are directly 
related to the radiative transfer properties of a compact Comptonizing region close to a 
black hole (e.g., Chakrabarti and Titarchuk 1995, CT95; Ling et al. 1997; Chakrabarti 
& Manickam, 2000, CM00; Muno et al., 2000; Feroci et al. 2000; Homan et al., 2001). 
These shocks have been studied most extensively in the literature (Chakrabarti 1989, C89; 
Chakrabarti 1990, C90) and their properties have been verified by several independent groups 
(Yang and Kafatos 1995; Nobuta and Hanawa 1994; Lu and Feng 1997) of workers. If for a 
given set of initial parameters the standing shocks exists analytically, numerical simulations 
also would find them (Chakrabarti and Molteni 1993; Molteni, Lanzafame and Chakrabarti 
1994; Ryu et al. 1995; Molteni, Ryu and Chakrabarti 1996) otherwise, the shock would be 
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oscillating (Molteni, Sponholz and Chakrabarti 1996; Ryu, Chakrabarti and Molteni 1997) 
causing QPOs in the X-rays emitted from the post-shock region. Furthermore, it has been 
observed that at least in some of the black hole candidates, such as GRS1915+105, outflows 
are produced from the region which also emits the Comptonized photons (Fender et al 2000; 
Dhawan et al. 2000) i.e., the post-shock region or the centrifugal barrier dominated region 
according to our present understanding of the flow solutions. 

It is therefore pertinent to ask if the properties of the shocks, such as location, strength, 
compression ratio etc. could be understood solely analytically. The solutions obtained so far 
by Chakrabarti and his collaborators and other groups always resort to numerical means. 
Roughly, the method was the following: for a given set of parameters (such as the specific 
energy and specific angular momentum), it is first determined if the flow allows more than one 
X-type sonic points (C89). This is because, at the horizon, matter must have velocity same as 
that of the velocity of light and for causality arguments, the flow must be supersonic. Thus, 
once the flow becomes supersonic through the outermost sonic point, and forms a shock 
(i.e., jumps to a subsonic branch) it must pass through the inner sonic point to become 
supersonic at the horizon. Next question is whether the specific entropy at the inner sonic 
point is higher compared to that at the outer sonic point. This is because, at the shock, 
entropy must be generated and post-shock flow which passes through the inner sonic must 
have higher entropy. The final and the most important question is whether all the Rankine- 
Hugoniot conditions (also known as the shock conditions) are satisfied somewhere between 
the two X-type sonic points. The three Rankine-Hugoniot conditions could be combined 
to obtain a combination of Mach numbers (Shock Invariant). C89 obtained this expression 
which is continuous across the shock and used this to obtain the shock location by iterative 
technique. 

In the present paper, we follow the same philosophy, but obtain the shock locations 
analytically. There were two motivations to do this. First, from the theoretical point of 
view, it is challenging to find solutions of a large number of non-linear equations which must 
satisfy a number of conditions mentioned above. Second, from the observers' point of view, 
any observation which could require standing shock waves, could be explained using more 
fundamental parameters, such as specific energy and angular momentum or even better, using 
accretion rates of Keplerian and sub-Keplerian flows as in a two-component flow solution of 
CT95. Analytical work also gives insight into why the shocks form in the first place. The 
boundary of the parameter space for shock formation is also obtained by analytical means. 
These findings are important as they would tell us when QPOs may or may not be seen. 
These details would be discussed in future. 

In the next Section, we present the model equations and shock conditions. In §3, we 
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present the sonic point analysis and write down the expression for Shock Invariant (C89) for 
the shake of completeness. In §4, we present the analytical expressions of the sonic points 
and discuss how the parameter space is divided into regions of one or three sonic points. 
In §5, we present expression of shock locations and again divide the parameter space into 
regions which may or may not have shocks. We compare our solutions with numerical work. 
In §6, we discuss some of the astrophysical implications of our solutions and finally, in §7, 
we make concluding remarks. 



2. Model Equations and Shock Conditions 

We start with a thin, axisymmetric, inviscid, steady flow which is in equilibrium in a 
direction transverse to the flow. The model equations which govern the motion of matter 
accreting into black hole are given as follows (C89): 

A. Radial Momentum Equation: 

dd_ \dP_ A 2 1 

dx p dx x 3 2{x — l) 2 

In a non- viscous flow, integration of this equation leads to the Energy Conservation Equation 

i9 2 n 2 \ 2 

^f + ^T + ^ + 9W (2) 

where, g(x) is the pseudo-Newtonian potential introduced by Paczyhski & Wiita (1980) 
which is given by, g(x) = — 2 <x-i) - Here, p is the mass density, P is the isotropic pressure, i? e 
and a e are the non-dimensional radial and the sound velocities measured in units of velocity 
of light c, x is the non-dimensional radial distance measured in units of the Schwarzschild 
radius r g = 2GM/c 2 , M being the mass of the black hole and G being the gravitational 
constant. 7 is the adiabatic index of the flow and P = Kp 1 is assumed to the equation of 
state. K is the measure of the specific entropy which is constant except at the shock location 
where local turbulence generate some entropy enabling the flow to pass through the inner 
sonic point. The subscript e refers to the quantities measured on the equatorial plane. 



B. Continuity Equation: 

■^-($ e pxh) = (3) 

ax 
which is integrated to obtain the Mass Conservation Equation : 

M = $ e pxh, (4) 
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where, h is the half-thickness of the flow at radial coordinate x. Here M is the mass accretion 
rate apart from a geometric constant. If we assume that the flow is in hydrostatic equilibrium 
in the transverse direction, the vertical component of gravitational force balances the pressure 
gradient force. Hence the expression for the half thickness of the disc is given by (C89), 

h(x) = a e x 1 ' 2 (x — 1). (5) 

We write the mass flux conservation equation in terms of i? e and the sound speed on the 
equatorial plane a e = \J^P/ p in the following way, 

M = $ e a q e x 3/2 (x - 1) = Aa q J(x) (4a) 

where, q = ^^ and f(x) = x 3 ' 2 (x — 1). 

The shock conditions which we employ here are the following (C89): 

(a) the energy flux is continuous across the shock — 

£+ = £-, (6) 

(b) the mass flux is continuous across the shock — 

M+ = M_ (7) 

and finally, (c) the momentum balance condition — 

W + + E+^ + = W _ + S _^2_ ( 8 ) 

where subscripts "— " and "+" refer, respectively, to quantities before and after the shock. 
A shock satisfying these conditions is termed as a Rankine-Hugoniot shock (Landau and 
Lifshitz 1959). Here, W and X denote the pressure and the density, integrated in the vertical 
direction [see, e.g., Matsumoto et al. (1984)], i.e., 

h 

S = / pdz = 2pJ n h, (9) 

-h 

and 

h 

W = f Pdz = 2P e I n+1 h, (10) 

-h 

where, /„ = L "'A, , n being the polytropic index as defined previously. In the subsequent 
analysis we drop the subscript e if no confusion arises in doing so. 
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3. Sonic Point Analysis and Shock Invariants 

In order to have a shock, the flow must be supersonic i.e., the stationary flow must pass 
through the sub-sonic flow to a supersonic flow. Discussions in this Section are based on 
early works [see, Chakrabarti (1989, 1990)]. 



3.1. Sonic Point Conditions 

From the radial momentum equation and the continuity equation (Eqs. 1 and 3) we 
derive the sonic point condition (or critical point condition) in the usual way (C89). The 
first derivative of the radial velocity w.r.t radial distance is given by, 

2a 2 dln(f) _ dG 

(11) 



u>" 7+1 dx dx 



dx 



•d 



2a 2 



Here, G(x) = -^ — 2 <x-i) * s ^ ne effective potential. Since the denominator must vanish at 
the sonic points, if the flow is assumed to be smooth everywhere, the numerator must also 
vanish simultaneously. The vanishing of the denominator gives, 

<?M = j^aKxJ. (12) 

The factor (which is unity only in isothermal flows) in front of a%(x c ) arises because the flow 
is assumed to be in vertical equilibrium. The vanishing of the numerator gives, 

,_ ( 7 + l)(x c -l)[A^(x c )-A 2 ] 
aAXc) ~ x\ (5z c -3) • [ ' 

The subscript c denotes quantities at the critical points. Here Xk is the Keplerian angular 
momentum defined as X K = x\j\l(x c — l) 2 ]. It is to be noted that since the left hand 
side is always positive, angular momentum at the sonic point must be sub-Keplerian, i.e., 
X(x c ) < Xk (e.g., C90 and references therein). 



3.2. Mach Number Relation at the Shock 

From the equations given in §2, we now seek an invariant relation which must be satisfied 
at the shock (C89). We rewrite the condition of energy flux continuity (Eq. 6), and the 
pressure balance condition (Eq. 8) in terms of the Mach number M = d/a of the flow, 

\M\a\ + -^ = \mW. + ^ I , (14) 
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M + = M + aif(x s ), (15) 

M- = M_aif{x s ), (16) 



where v = — ^r, arid 



+ Mi ) = -^— + M 1 , (17) 

where, v = -^- and x s is the location of the shock. From Eqs. (14-17) one obtains the 
following equation relating the pre- and post-shock Mach numbers of the flow at the shock 
(C89), 

_ [M + (3 7 - 1) + (2/M + )] 2 _ [M_(3 7 - 1) + (2/M_)] 2 

2 + ( 7 -l)M 2 2 + ( 7 -l)M£ ' { } 

The constant C is invariant across the shock. The Mach number of the flow just before and 
after the shock can be written down in terms of C as, 



H3r,-l)-C±y/C^ m (19) 

(7-l)C-(3 7 -l) 2 ' UH ' 

The product of the Mach number is given by, 

M + M_ = m . (20) 

[(3 7 -l) 2 -( 7 -l)C] 1/2 

4. Analytical Expression of the Sonic Points and Behavior in Parameter Space 

To obtain shock locations, we need to first obtain the locations of the sonic points, and 
ensure that at least two of them are X-type (C90). In §3.1, we presented the sonic point 
conditions. Using the definition M = ft/a of the Mach number, and substituting a^(x c ) from 
Eq. (13) into Eq. (2), we obtain following algebraic equation for x c given by, 

Mx A c - Ox\ + Vx\ -Qx c + TZ = 0. (21) 

where, M = 105; O = 16£ + 2n - 3; V = 6£ + A 2 (4n - 1) - 3; Q = 8n\ 2 ; K = (1 + 4n)A 2 
and n (= -^-) is the polytropic index. 

We solve this equation analytically (Abramowitz and Stegun 1970) and get the location 
of the sonic points. Details are given in Appendix A. For the purpose of critical or sonic 
points, D of the Appendix would be denoted as D c . The equation has four roots and Eq. 
A. 8 can be used to check whether all of them are complex (pair of complex conjugates) or at 
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least two are real and two are complex or all four are real. A necessary condition to form a 
shock wave is to have all four real roots. Of course, only one would be inside the black hole, 
and the other three would be outside, and out of these, due to topological reasons, only two 
would be X-type or saddle type and the one in between must be 'O'-type or Center type. 
These are determined by computing derivative ddjdx at the sonic point by using L 'Hospital's 
rule and checking if they are real. For our purpose, two derivatives at each sonic point must 
be real and of opposite signs in order that the sonic point be of X-type or saddle type. 

Figure 1 shows the division of the parameter space. Denoting the discriminant D (of 
Appendix A) by D c we find that the condition D c < is the necessary condition to have 
three sonic points. The boundary D c = separates this region on which two sonic points 
merge and the third one remains separate. Outside of this region D c > 0, and only one sonic 
point is possible and the other two roots are complex conjugate of each other. The dotted 
curve in the middle represents the condition: 

Mo = M h 

where Aii and Ai are the entropy accretion rates at the inner and outer sonic points 
respectively. The region above it contains parameters with Aii < M. Q and the region below 
it contains parameters with Aii > hA . This latter region is suitable for shock formation in 
accretion flows. 

In passing we wish to point out that from Eq. (11) one can easily show that locations 
where dM/dx = exactly coincide with the sonic points of the flow. Thus the number of 
extrema of M — M(x) is the same as that of the sonic points and hence the division in Fig. 
1 could give an idea about the behavior of M — M(x) as well. 



5. Analytical Expression for the Shock Location and Behavior in Parameter 

Space 

A black hole accretion flow being transonic, it must satisfy two sonic point conditions at 
the cost of one extra unknown, namely, the the sonic point. Because of this extra condition, 
out of the three constants of motion, namely, £ , M and A, only two are to be supplied as 
free parameters. While computing shock locations, we used past experience derived from 
numerical methods, that only one of the shocks, namely x s ^ (C89 notation) is stable, and 
accordingly our procedure as delineated below attempts to compute only this location. 

The flow will have a shock only at the point where the shock invariant condition is 



satisfied. Simplifying the shock invariant relation (Eq. 18), we obtain, 

2( 7 - l)(Ml + M 2 _) - [(3 7 - l) 2 - 2(3 7 - 1)( 7 - 1)]M 2 M! +4 = 0. (22) 

We consider a relativistic flow with 7 = 4/3 so that the polytropic index n — 3. Then from 
Eq. (22), we obtain, 

2(M 2 + M 2 ) - 21M 2 M 2 + 12 = 0. (23) 

We now expand the post-shock Mach number M 2 by a polynomial which must satisfy 
the following conditions: 

(a) ^p is zero at the central 'O'-type sonic point. This is a general property of the flow 
(vide discussion at the end of last Section). 

(b) Mach no. (M + ) at the location of the middle sonic point must match with that 
derived from approximate analytical solution obtained using the energy equation (Eq. 2 and 
4a). 

(c) Solution must pass through the position where ^ is 00 and a good guess of this 
the location (say, from the location of the sonic points) is known. 

Similarly, we expand M 2 by a polynomial which must satisfy the following conditions: 

(a) ^f is zero at the outer sonic point (location of which is already determined above). 
This is a general condition (vide discussion at the end of the previous Section). 

(b) Mach No. at the outer sonic point must match with the analytical value obtained 
from the Sonic point condition (Eq. 12). 

(c) Mach No. (M_) at the location of the middle sonic point must match with the 
analytical value derived from approximate analytical solution obtained using the energy 
equation (Eq. 2 and 4a). 

Keeping in mind that an algebraic equation which is beyond quartic cannot be solved 
analytically [see, Abramowitz and Stegun (1970)], we expand M\ as quadratic equations so 
that Eq. (23) may become quartic. We shall show a posteriori that such an assumption 
introduces a very small and tolerable error in our computation. 



If x s denotes the shock location, we assume, 

2 
M± = £4*±F2 (24) 



where, ^4[<j,±i are constant co-efficients to be determined from the conditions mentioned above. 
We find them to be, 
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i-(mI) 

V ' / mid 
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[2,+]> 



\Mi) -1 

V / mid 

Kp^out *&mid) 

—lX, m i(iA\2~ 
lut A [2,- 



1 + X nn, + A 



where, x mi d, x in f and x out are the middle (O-type) sonic point, the position where first 
derivative of Mach No. is infinity and the outer sonic point respectively. 

We now substitute the above expression (Eq. 24) for the Mach number in the Mach 
invariant relation (Eq. 23) to obtain the following algebraic equation: 



<J~\, Uj „ | & dj „ | L^ dj , 



Vx s + T = 0, 



(25) 



where, 

y = (3 7 -l)( 7 + l), 
Z = 2(7-1), 
A = yA [2j+] A [2 _ ] , 
B = y (A [li+] A [2iH + 
c = y (A [0t+] A [2H 
v = y (A [0:+] A [1: _ 



A [2,+] A [1,-]) > 

Al,+] Al "I + A 2 ,+]A° ,-]) - -2 (-^[2,+] 

V]A[o r ])-2(i[i, + ]+i[i rl ), 



^p,-])» 



T = yA [QM A [Q ^ - Z (A [0i+] + A [0 ,_]) - 4. 

We solve this equation for X s analytically using the same procedure as in §4 (details 
in Appendix A). We denote the discriminant D by D s , and Q and R values as Q s and R s 
respectively for our discussion on the parameter space behavior of shocks. 

In Figure 2, we redraw the parameter space as in Figure 1, but consider the formation 
of shocks alone. We find that Q s > produces no shock from above, and Q s = with 
R s 7^ gives the boundary of the weakest shock (shocks with unit compression ratio). This 
boundary, though obtained using our approximate analytical method, generally coincides 
with the dotted curve of Figure 1. The edge of the boundary is obtained with an extra 
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condition R s = 0. Thus, R s progressively decreases towards the edge along the dashed 
curve. This edge (denoted by D s = R s = 0) ought to have coincided with the cusp of 
the D c = curve drawn for the sonic point (see, Fig. 1 also), had the analytical method 
been exact. A small shift is the evidence that a small error is present at this corner of the 
parameter space. We also provide the region of the oscillating shocks (Q s < and D s < 0). 
Here, the shock location is imaginary and therefore shock continuously oscillates back and 
forth causing a very interesting astrophysical effect known as Quasi-Periodic Oscillations 
(QPOs). This would be discussed in the next Section. The boundary between the shock and 
no-shock region from below is denoted by the dashed curve marked by D s = 0. Below the 
no-shock region where the energy and angular momentum of the flow are very low, the flow 
has only one sonic point. 

In Figure 3, we draw shock location x s (along Y-axis) as a function of energy (along 
X-axis). Different curves are drawn for different specific angular momentum of the flow. 
The rightmost one is for A = 1.51 and the leftmost one is for A = 1.84, interval of A is 
0.01. As angular momentum is increased, the shock is located farther from the black hole. 
Comparing the locations from those obtained analytically, one notices that the same location 
is obtained for a specific angular momentum slightly more (~ 3%) than that used in the 
numerical method. We therefore believe that the results obtained are very much reliable. 

In Figure 4, we present a comparison of the boundary of the parameter space for which 
standing shocks may form as obtained by our analytical solution (shaded region) and by the 
numerical means (solid curve) existing in the literature (C89). The agreement is very good 
except in a region near the cusp (as also noted while discussing Fig. 2). Since very little 
parameter space is involved at this edge, we think that this small mismatch is tolerable. 



6. Astrophysical Applications 

Even though a black hole has no hard surface, it is remarkable that matter forms 
standing shocks around it in the same way as a shock is formed when a supersonic flow 
encounters a hard boundary. Shock waves heat up a gas and puff it up. This post-shock 
region intercepts soft photons from the pre-shock matter, particularly from the Keplerian 
disk which is located on the equatorial plane (CT95). In this scenario, the nature of the 
Comptonized radiation depends on the amount of matter in the sub-Keplerian and in the 
Keplerian flow: if the intensity of soft photon is very low, they cannot cool down the post- 
shock region by inverse Compton process and the spectrum remains very hot. On the 
other hand, if the intensity of soft photon is very high (i.e., the Keplerian rate is large), 
they cool down the post-shock region to the extent that the shock could not be sustained 
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(pressure balance condition breaks down). This produces a soft state spectrum with a hard 
tail due to bulk motion Comptonization (CT95). There are several models in the literature 
which perhaps explain the soft/hard states. However, no model other than CT95 explains 
the power-law hard-tail in the soft state. Similarly, when the question of Quasi-Periodic 
Oscillations (QPOs) come, shock-oscillation model turns out to be a sufficiently satisfactory 
one (CMOO). 

When the parameters fall in the 'no shock' region of Fig. 2, shock location becomes 
imaginary. However, three sonic points are still present and the entropy of the flow at the 
inner sonic point continues to be higher compared to that at the outer sonic point. In this 
case, shock starts oscillating with a time period (T s ) comparable to the infall time from the 
post-shock region (Ryu, Chakrabarti and Molteni 1997). Even when shocks form, if the infall 
time-scale turns out to be comparable to the cooling time, then the resonance condition is 
satisfied (Molteni, Sponholz and Chakrabarti 1996; Chakrabarti and Manickam 2000) and 
shocks oscillate in time scales of: 

T s ~ x s /v s , 

where, v s ~ R~ 1 x s is the infall velocity and R is the compression ratio at the shock (easily 
obtained analytically from our equations). Observed QPO frequencies are comparable to 
±/T s . 

When a mixture of Keplerian and sub-Keplerian matter is accreted, it is easy to obtain 
the parameters A and £ in terms of the Keplerian (M d ) and sub-Keplerian (Mh) accretion 
rates. Suppose, the viscosity parameter is such that the flow is deviated from a Keplerian disk 
at x — xk where its energy and angular momentum were £ d and A^ respectively. Further 
suppose that the sub-Keplerian halo has a constant energy £ ~ ~ £ h and constant angular 
momentum A/j, then the average angular momentum and energy of the transonic flow would 
be, 

M d £ d + M h £ h 



<£>- 
<X>-- 



M d + M h 
M d \ d + M h X h 



M d + M h 

It is easy to compute the shock location of the resultant flow using our formalism given 
above. 

It is to be noted from Fig. 3 that shock solutions are allowed only if the specific energy is 
positive. In other words, if a flow deviates from a cool Keplerian disk on the equatorial plane, 
the flow cannot have shocks as the specific energy in such a flow would be negative unless 
this flow is mixed with a substantial amount of sub-Keplerian matter with a positive energy. 
Typically, M^ >> M d and even with a small energy, the specific energy of the mixture 
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becomes positive, giving rise to shocks and (unbound) winds. In the case when magnetic 
dissipation is present, flow energy could increase to a positive value and a solution with a 
shock would be allowed. The prospect of magnetic energy dissipation has been discussed 
by several workers in the literature (Shapiro 1973; Bisnovatyi-Kogan and Blinnikov 1976; 
Bisnovatyi-Kogan 1998). Briefly, since the magnetic field rises as B r oc r~ 2 and therefore 
magnetic pressure rises as P ma g oc r~ 4 while the gas pressure in the sub-Keplerian matter 
goes as P gas oc r~ 5 / 2 , magnetic field in excess of the equipartition value would escape from 
the disk buoyantly and may dissipate at the atmosphere as in the case of the Sun. If the 
flow has specific energy Sh at, say r = 100r 9 where the flow was in equipartition, then at 
the shock, the energy would be at least 2 4 = 16 times larger if all the magnetic energy is 
dissipated into the flow. Thus, a basically free-fall matter of £ ~ 10~ 4 would have an energy 
~ 10~ 3 and thus a shock at a few tens of Schwarzschild radii would be expected. 

An accreting flow can intercept hard X-rays emitted at the inner edge and this pre- 
heating effect need not be negligible. For instance, a flow emitting isotropically with 6% 
efficiency will definitely intercept 9 ~ Gj„/47r fraction of radiation in between the shocked 
region and the Keplerian disk. Assuming 9 ~ 0.1, the energy deposition due to pre- heating 
is 0.006 which is significant. This would energize Keplerian matter as well and shocks in the 
sub-Keplerian flow would be expected. 



7. Concluding Remarks 

So far in the astrophysical literature, existence of shocks in accretion flows has been 
indicated by steady and time-dependent numerical simulations. Study of these standing and 
oscillatory shocks in an accretion flows has been shown to be of great importance. Presently 
we show that the shocks could be studied completely analytically, at least in the case of 
thin, axisymmetric, inviscid flows with positive energy. We note that shock locations vary 
with flow parameters in a simple way — they form farther from a black hole when angular 
momentum is increased. This proves that they are mainly centrifugal pressure supported. 

Given that the shocks, especially the standing shocks are the ideal locations at which a 
flow is heated up, hard X-rays are produced from the post-shock region after the soft-photons 
are processed by the flow by inverse-Comptonization process. Thus, the spectral states and 
the time-dependent behavior of the hard X-rays are directly related to the behavior of this 
region. For instance, CT95 computes steady state spectra with post-shock region as the 
source of hot, Comptonizing electrons. CM00 establishes that QPOs are due to oscillations 
of this region, since only hard X-rays are seen to exhibit QPOs. We therefore believe that 
shocks should be an important ingredient of an accreting system. However, if the disk is cool 
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and Keplerian very far away, the specific energy must be negative. So the problem is not 
whether shocks should exist, the problem seems to be how to energize Keplerian matter as 
it becomes sub-Keplerian, by, say, magnetic energy dissipation, pre-heating etc. These work 
is in progress and would be reported elsewhere. 

From the observers point of view, our work could also be useful, since the steady spectra, 
QPO frequencies, etc. are, in principle, determined analytically from a few free physical 
parameters. In future, we shall apply our understanding to obtain spectral properties and 
QPO behaviors more quantitatively. 

This work is partly supported by a project (Grant No. SP/S2/K-14/98) funded by 
Department of Science and Technology (DST). 
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A. Method of the Analytical Solution 

The procedure of obtaining an analytical solution of a quartic equation, 

q 4 + hq 3 + b 2 q 2 + b 3 q + b 4 = Al 

is to first obtain a solution of the following cubic equation: 

p 3 + axp 2 + a 2 p + a 3 = 0, A. 2 

where, 

«i = — b 2 , a 2 = &1&3 — 46 4 , and a 3 = 46 2 &4 — b\ — b\b 4 . 

This preprint was prepared with the AAS IATgX macros v5.0. 
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Let, 



p. 3a2— a\ p 9aia2—27a:j—2a'f 

V — 9 j tt — 54 ; 



S = y R + y/g 3 + R 2 , T=^R- VQ 3 + R 2 . 
The discriminant is defined as 



D = Q 3 + R 2 . A3 

If D > 0, one root is real and two roots are complex conjugate. In this case, the real solution 
is, 

Pi = 5 , + T--oi. AA 

o 

If D — 0, all roots are real and at least two are equal. 
If D < 0, all roots are real and unequal. They are: 

l„s 1 



P2 



pi = 2y/-Qcos(-9) - -d, A5 



2 v /^gcos(^ + 120°)-ia 1 , A6 



and, 



where, cos^ = R/^—Q. 



P3 = 2 v /Z Qcos(^ + 240°)-^a 1 , A7 



One can now write a quadratic equation using any one of the real solutions of the cubic 
equation (see, Spiegel, 1968) as follows: 



.2 1 



z 2 + -{h ± yJVi - 46 2 + 4 Pl }z + -{ Pl T \Jvl - 46 4 } = 0. A8 

This is a quadratic equation which can be solved easily. Since we applied this procedure 
both of the sonic points and shocks, we denoted quantities like, D,Q and R by D C ,Q C and 
R c for sonic (critical) points and D S ,Q S and R s for shocks respectively. 
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Fig. 1. — Division of the parameter space as spanned by the pair £ , A according to the number 
of sonic points. Solid curve represents D c = which divides the region into D c > (one 
sonic point) and D c < (three sonic points) regions. The plot is for 7 = 4/3. For 7 £ 1.5, 
D c > always, suggesting that no shocks are possible in a flow with vertical equilibrium. 
The dashed curve further divides the region into two regions where entropy accretion rate 
Ai at the two saddle type sonic points behave differently (inner point is denoted by i and 
outer point is denoted by o). 

Fig. 2. — Division of the parameter space as spanned by the pair S, A according to whether 
shocks would form or not. Solid curve represents D c = as in Fig. 1. Dashed curve (D s = 0) 
surrounds the region with shocks in accretion. When D s < and yet, there are three sonic 
points, shocks are oscillatory, giving rise to quasi-periodically varying hard X-rays. 

Fig. 3. — Variation of shock location (x s along y-axis) with specific energy (S along x-axis) of 
the flow. Each curve is drawn for a specific angular momentum A. From right to left curved 
are drawn for A = 1.51, 1.52, 1.53,... till 1.84 respectively. For a given specific energy S, 
shock location increases with increasing centrifugal force (through A). Similarly, for a given 
A, shock location increases with energy. 

Fig. 4. — Comparison of the boundary of the parameter space in the (£, A) plane using 
the numerical and the analytical methods. Except for the region near the cusp (upper left 
corner) the agreement is very strong. 
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